CHARACTERISTIC MATRICES FOR LINEAR PERIODIC DELAY DIFFERENTIAL 

EQUATIONS 

JAN SIEBER* AND ROBERT SZALAI+ 

Abstract. Szalai et al. (SLAM J. on Sci. Comp. 28(4), 2006) gave a general construction for characteristic 
matrices for systems of linear delay-differential equations with periodic coefficients. First, we show that 
matrices constructed in this way can have a discrete set of poles in the complex plane, which may possibly 
obstruct their use when determining the stability of the linear system. Then we modify and generalize the 
original construction such that the poles get pushed into a small neighborhood of the origin of the complex 

T— H plane. 
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^~~^ 1. Introduction. Linear delay differential equations (DDEs) with coefficients pe- 

^— I riodic in time come up naturally when one analyses nonlinear DDEs: one of the most 

^N common and simplest invariant sets encountered in nonlinear systems are periodic 

^ orbits [5] . If one wants to find out if a periodic orbit is dynamically stable or unstable, 

C/^ or how its stability changes as one changes system parameters one has to linearize 

Q the nonlinear DDE system in the periodic orbit, which results in a linear DDE with 

• time-periodic coefficients [2, 6] . This linear DDE is typically written in the form 

d x(,t) = L(,t)x, (1.1) 

s 

i_i where the coefficients (hidden in the operator L^t)) are periodic in t. Without loss of 

generality we can assume that the time dependence of L is of period 1 (this corresponds 
\^ to rescaling time and L). Then the exponential asymptotic stability of the periodic 

p^ orbit in the full nonlinear system is determined by the eigenvalues of the time-1 map 

psj T generated by the linear DDE (1.1). See the textbooks [2, 6] for a summary of 

lO the fundamental results on Floquet theory and [11] for robust numerical methods 

'^ to calculate these eigenvalues in practical problems. For analytical purposes (and, 

l/~j possibly, for numerical purposes) it is useful to reduce the eigenvalue problem for the 

^D infinite-dimensional time-1 map T to an algebraic problem via a characteristic matrix 

^^ A(A) e C"^". One would expect that this matrix A should satisfy, for example, that 

. . (i) A(A) is regular if (and only if) A e C is in the resolvent set of T (that is, 

^ A/ — T is an isomorphism), 

(ii) A is a root of det A(-) of order q if and only if A is an eigenvalue of T of 
algebraic multiplicity q (see [8] or Lemma 4.2 in §4 how one can construct the Jordan 
Cd chains of T from A), and 

(iii) dimker A(A) is the geometric multiplicity of A as an eigenvalue of T. 
An immediate consequence of the existence of such a matrix A would be the upper limit 
n (the dimension of A(A)) on the geometric multiplicity of eigenvalues of T. 

For time-invariant linear DDEs (where L does not depend on t in (1.1)) the exis- 
tence of a characteristic matrix has been known and used for a long time. The general 
theoretical background and construction (extended to more general evolutionary sys- 
tems such as, for example, neutral equations and some classes of hyperbolic partial 
differential equations) is given by Kaashoek and Verduyn-Lunel in [8]. 

For time-periodic linear DDEs the textbook of Hale and Verduyn-Lunel [6] has 
developed characteristic matrices only for the case where all time delays are multiples 
of the period. This has been generalized to delays depending rationally on the period in 
order to derive analytical stability results for periodic orbits of the classical scalar delayed 
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positive feedback problem [12]. A characteristic function in the form of a continued 
fraction expansion is given for a scalar DDE in [7]. In general, precise analytical 
statements about the location of Floquet multipliers (or the number of unstable Floquet 
multipliers) can be made only for DDEs with very special structure, typically scalar DDEs 
with single delays and particular requirements on the right-hand side, see [1, 3, 14, 15] 
for classical works. An extension of these results to cyclic DDEs using a discrete Lyapunov 
function is given in [10]. 

Szalai et al. give a general construction for a characteristic matrix A for time- 
periodic DDEs in [13] that lends itself easily to robust numerical computation. Since 
one expects that A(A) must have an essential singularity at A = (the eigenvalues of 
the compact time-1 map T accumulate at 0) the constructed matrix is rather of the 
form A(/x) where /i = A~^. We show (in §3) that the characteristic matrix A(/^) as 
constructed by Szalai et al. typically has poles in the complex plane, making it unusable 
whenever eigenvalues of T of interest coincide with these poles. However, the set 
of these poles is discrete and accumulates only at oo. Moreover, in §4 we provide a 
modification A;t(/i) of the original matrix A(/>(), which allows us to push the poles of 
A;^ to the outside of any given ball of radius R > 1 in the complex plane by increasing 
the dimension of Aj^. This modification reduces the eigenvalue problem for T to a 
root-finding problem for det A^(/i) for all eigenvalues of T with modulus larger than 
1/R, which is useful because one is typically interested in the largest eigenvalues of T. 
To keep the notational overhead limited we develop our arguments for the case of a 
DDE with a single delayed argument and a constant delay t < 1. In Appendix B we 
show that the generalization to arbitrary delays (distributed and reaching arbitrarily far 
into the past) is straightforward. 

One immediate application for the characteristic matrix came up in [16] where 
Yanchuk and Perlikowski consider periodic orbits of nonlinear DDEs with a single fixed 
delay t but change the delay to t + fcP (where k is an integer and P is the period 
of the periodic orbit), and then study the stability of the periodic orbit (which does 
not change its shape) in the limit fc ^ oo. The limit is, of course, a singular limit if 
one considers the time-P map T but the characteristic matrix A^^ has a well-defined 
regular limit. This permitted the authors to draw conclusions about the linear stability 
properties of periodic orbits for sufficiently large delays from the properties of the 
so-called pseudo-continuous spectrum [16]. As [16] relies on the original construction 
from [13], which may have poles near the unit circle, our note closes a gap in the 
argument of [16]. 

2. Construction of characteristic matrix for a single fixed delay. Consider a 
periodic linear differential equation of dimension n with a single delay t and continuous 
periodic coefficient matrices j4(t) and B[t) (both of size nx n): 

x(t)=A(t)x(t)-|-B(t)x(t-T) (2.1) 

where we assume that the period of the time dependence of ^4 and B equals 1 without 
loss of generality (this corresponds to rescaling time, t, A and B). We also assume for 
simplicity of notation that t < 1. Let us denote the monodromy operator (also called 
the time-1 map) for (2.1) by T : C([-1,0];C") -^ C([-1,0];C"). C([-1,0];C") is the 
space of continuous functions on the interval [—1,0] with values in C". That is, T maps 
an initial condition x in C([— 1,0];C") to the solution of (2.1) after time 1 starting 
from the initial value x(0) (the head point) and using the history segment x(-). For 
X e C([-1,0];C") we define Tx precisely as the solution y e C([-1,0];C") of 

y(0=^(tMO + B(t)^J^-^) , ift.[.-l,0] 

\x(l-Ft-T) ifte[-l,T-l) (2.2) 

y(-l) = x(0). 
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The complex number A is an eigenvalue of T (a Floquet multiplier for (2.1)) if one can 
find a non-zero x such that Tx = Pix. We know that T is compact [6] for t < 1 (and 
that, for arbitrary t, we can find a power T™ that is compact). Hence, any non-zero A 
is either in the resolvent set of T (that is, XI — T is regular) or it is an eigenvalue of T 
with a finite algebraic multiplicity. 

If we introduce p, = 1/A we may write the eigenvalue problem for T as follows: 

.... r.ix[t-T) ifte[T-l,0], 

x[t)=AMx[t) + B[t)\ \ J ' '' (2.3) 

[fxx[l + t-z) if t e [-1,T- 1), 

x(-l) = /ix(0) (2.4) 

where x e C([— 1,0];C") is continuous. Reference [13] constructs a characteristic 
matrix A(/i) for T based on the following hypothesis. 

Hypothesis 2.1 (Szalai'06) 

The initial- value problem on the interval [—1,0] 

.... ..f^(f-T) ifte[T-l,0], 

x(t)=A(Ox(t) + B(0 \ J ' ^' (2.5) 

[lj,x[l + t-T) if t e [-1,T- 1), 

x(-l) = v (2.6) 

has a unique solution xeC([-l,0];C") for all /xeC and v eC". 

Note that the problem (2.5)-(2.6) differs from the boundary-value problem (2.3)-(2.4) 
because we have replaced the boundary condition x(— 1) = fxx[0) by an initial condition 
x(— 1) = V. The formulation of the hypothesis in [13] is (only apparently) weaker: it 
claims the existence only for /x e C for which l//i is either in the resolvent set of T or 
an eigenvalue of T of algebraic multiplicity one. If Hypothesis 2.1 is true then one can 
define the characteristic matrix A(/i) via 

A(/i)v = V - /ix(0) (2.7) 

where x(0) is the value of the unique solution of (2.5)-(2.6) at the end of the time 
interval [-1,0]. 

Integral equation for (2.5)-(2.6). For system (2.5)-(2.6) one has to clarify what it 
means for x e C([-1,0];C") to be a solution. We call x e C([-1,0];C") a solution of 
(2.5)-(2.6) if X satisfies the integral equation 

x(t) = v-l-Mi(/x)[x](t) where (2.8) 

[Mi(M)x](0= f'A(s)x(5) + B(5)^^';^^ ^ '!'^!^7^'°l\d^ (2.9) 

J_j [jL6x(l-|-S - Tj lfse[-l,T-lJ 

pointwise for all times t on the interval [—1,0]. The linear operator M^C/i) maps 
C([-1,0];C") back into C([-1,0];C") and is compact. (Mj will be generalized later 
toM^ forfc > 1.) 

Lemma 2.2 

For fi sufficiently close to the initial-value problem (2.5)-(2.6) has a unique 
solution for all v. 

Proof:. For (2.8) to have a unique solution it is enough to show that / — M^du) is 
invertible. Since the time-1 map T is well defined and linear in x = (that is, Tx = 
for X = 0) we have that (2.2) has only the trivial solution for x = 0. This in turn implies 
that the corresponding integral equation 

y = Mi(0)y (2.10) 
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has only the trivial solution in C([— 1,0];C"). Thus, the kernel of/ — Mi(0) is trivial, 
and, since / — Mi(0) is a compact perturbation of the identity, / — Mi(0) is invertible. 
Moreover, the real-valued function /i -^ ||Mi(/^) — Mi(0)|| (using the operator norm 
induced by the maximum norm of C([— 1,0];C")) is continuous in n because Mj 
depends continuously on ji (see definition of M^ in (2.9)). This implies that / — M^C/i) 
is also invertible for small /i, which in turn means that (2.8)-(2.9) has a unique solution 
for small /i. D 

The following lemma states that the function on the complex domain /i -^ [/ — 
Mi(/i)]"\ which has operators in if(C([-l,0];C");C([-l,0];C")) as its values, is 
well defined and analytic everywhere in the complex plane except, possibly, in a discrete 
set of poles. 

Lemma 2.3 

The operator-valued complex function fx-^ [I — Mi(/x)]~^ can have at most a 
finite number of poles in any bounded subset of C . All poles of [/ — Mj (/i)] ~ ^ (if 
there are any) are of finite multiplicity. If /x is not a pole then /x -^ [/ — Mi(/>i)]~^ 
is analytic in /x. 

Proof:. We can split the operator MjC/i) into a sum of two operators: 

Mi(/i) = Mi(0)-l-|ai- where (2.11) 

[L.](0=f^w{° ^ , ';^\'-'''\6s. (2.12) 

J_i (x(I-I-s-t) ifse[-l,T-l) 

The operator / — Mi(0) is invertible (and a compact perturbation of the identity) and 
L : C([-1,0];C") ^ C([-1,0];C") is compact. Thus, we can rewrite / - MjC/x) as 

/ - Ml (/x) = (/ - Ml (0)) [/ - /x(/ - Ml (0))- 1 L ] and, thus, 

[/ - Mi(/x)]-i = A • [A/ - (/ - Mi(0))-1l] "' (/ - Mi(0))-i 

(keeping in mind that /x = 1/A). The function A -^ [A/ - (/ - Mi(0))"^L]"^ on the 
right-hand side is the standard resolvent of the compact operator (/ — Mi(0))~^L, which 
has only finitely many poles of finite multiplicity outside of any neighborhood of the 
complex origin according to [9], and is analytic everywhere else. D 

Lemma 2.3 carries over to the characteristic matrix A(/x) defined in (2.7): since 
X = [I — Mi(/x)]~^v (if we extend v e C" to a constant function), A(/x)v = v — /ix(0) is 
an analytic function with at most a discrete set of poles of finite multiplicity (possibly 
accumulating at oo) . 

The characteristic matrix A(/x) (where it is well defined) has the properties one 
expects: for a complex number /x in the domain of definition of A, l//x is a Floquet 
multiplier of the DDE (2.1) if and only if A(/i) is singular (det A(/i) = 0) . Any vector v in 
its kernel is the value of an eigenfunction for l//x at time t = — 1. The full eigenfunction 
can be obtained as the solution of the initial- value problem (2.5)-(2.6). If A(|U) is 
regular then l//x is in the resolvent set of the eigenvalue problem (2.3)-(2.4). The 
Jordan chains associated to eigenvalues 1/n can be obtained by following the general 
theory described in [8]. We give a precise statement (see Lemma 4.2 in Section 4) about 
the Jordan chain structure after discussing H5Apothesis 2.1. 

The construction can be extended in a straightforward manner to multiple fixed 
discrete delays and distributed delays. The only modification is that higher powers 
of /x have to be included for delays larger than 1. For example, if t e (1, 2] then the 
initial-value problem (2.5)-(2.6) is modified to 

.^ .^ , , f/x^(l + f-T) ifte[T-2,0], 

x(t)=A(0x(0+B(0]'^2 \ -f -, o, 

(/x^x(2-l- 1- t) if t e [-1,t -2), 
x(-l) = v. 
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Thus, also for arbitrary delays A(|U) always has at worst finitely many poles of finite 
multiplicity in any bounded subset of C. Appendix B discusses the case where the 
functional accessing the history of x is a Lebesgue-Stieltjes integral. 

3. The poles of the characteristic matrix. Hypothesis 2.1 is not true in its stated 
form (and neither in the form stated in [13]). Let us analyse the simple scalar example 

x[t) = x(t-^^ (3.1) 

on the interval [—1,0]. Problem (2.5)-(2.6) restated for this example is equivalent to 

Xi(t) = /XX2(t), Xi(0) = V 

X2(t) = Xi(t), X2(0) = Xi(l/2) 

where Xi and Xj are in C([0, 1/2];C), and Xi(t) = x(t - 1) and X2(t) = x(t - 1/2) for 
t in the interval [0, 1/2]. The time-1/2 map for the IVP of the ODE (3.2) with initial 
conditions x^CO) = v, X2(0) = w is 



Xi(l/2)" 

X2(l/2) 



coshQy/^) y/IsinhQy/l) 
sinhQy]u)/v7I cosh(iy]u) 



where -//x refers to the principal branch of the complex square root (the set of solutions 
is the same for both branches as all entries in the matrix are even). The continuity 
condition X2(0) = Xi(l/2) in (3.2) is satisfied if and only if 



1 — y/Isinh I — v7^ 1 



w + cosh ( - y/I I V = 0, 



which has a unique solution for arbitrary v if and only if 

Ot^ 1 — v^sinh j -^/Jji I . 

Thus, we see that for example (3.1) the initial-value problem (2.5)-(2.6) cannot be 
solved for non-zero v whenever fxisa root of the function 

/(m)= 1- vT^sinhf 2^^) • 

This function / is a globally defined real analytic function of /x (the expression on the 
right-hand side is even in yp) . It has an infinite number of complex roots accumulating 
at infinity. One of the roots is real: /(O) = 1, and /(/x) -^ — oo for /x -^ -l-oo and /x e K 
(let's call it ^q: /Xq ^ 1.8535). Consequently, example (3.1) provides a counterexample 
to Hypothesis 2.1. If one extends example (3.1) by a decoupled equation that has its 
Floquet multiplier at l//Xo (for example, y = — [log/ioJy) then it also contradicts the 
hypothesis as stated in [13]. 

As we have discussed in Section 2, one finds that /x = 1/A is a pole of A(/i) if A is 
a non-zero eigenvalue of the operator [/ — Mi(0)]~^L where M^ and L are defined in 
(2.9) and (2.12). Thus, whenever [/ — M^(0)]~^L has a non-zero spectral radius we 
must expect that A(/x) has poles. 

4. The extended characteristic matrix. The construction of A(/x) suffers from 
the problem that the poles of A(/x) may coincide with the inverse of Floquet multipliers 
of (2.1) of interest. A simple extension of the construction of A(/x) permits one to push 
the poles to the outside of a circle of any desired radius R. Then the characteristic matrix 
can be used to find all Floquet multipliers outside of the ball of radius 1/R. We explain 
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the extension in detail for a single delay t < 1 (the straightforward generalization to 
the general case is relegated to Appendix B) . 

The idea is based on the observation that the linear part of the right-hand side 
in the integral equation formulation (2.8)-(2.9) of initial-value problem (2.5)-(2.6) 
has a norm less than 1 if the interval is sufficiently short (thus, making the fixed point 
problem (2.8)-(2.9) uniquely solvable). 

Similar to a multiple shooting approach, we partition the interval [ — 1,0] into k 
intervals of size 1/fc: 



Ji= [ti,ti+i) 



i i + l\ 

-1-1-7,-1-1-—— fori = 0,...,fc-l. 
k k 



(4.1) 



Then we write a system of k coupled initial value problems for a vector of k initial values 



x(t)=A(t)x(t)-|-B(t) 



x(t-T) 



if te [t-1,0], 



^fjLx[l + t-T) if t e [-1,T- 1), 
x(t;) = V; for i = 0, ...,fc — 1 



(4.2) 
(4.3) 



where t e [—1,0]. System (4.2)-(4.3) corresponds to a coupled system of fc differential 
equations for the fc functions x\j. in the time intervals J;. More precisely, x is defined as 
a solution of the fixed point problem 



t) 



ifse [t-1,0]. 



x(t) = v,(t)+ I A(s)x(s) + B(s)r^^ "^ 7 ;- "~;'ds where (4.4) 

[lJLx{l+s — T) ifse[— 1,T — IJ 

v,(t) = V, if t e J„ (4.5) 

a^[t) = ti ifteJ;. (4.6) 



We point out that a solution x of (4.2)-(4.3) is not necessarily in C([-1,0];C") 

2.-(t2)- 




Fig. 4.1: Illustration of solution x e Q for multiple-shooting problem [4.2)-[4.3) with fc = 3 
sub-intervals. The right-sided limits x(t,)+ are equal to the variables V;. The gaps between the 
left-sided limits x(t,)_ and x(t,)+ are zero if [vq, ...,v,j] is in the kernel of the characteristic 
matrix Aj.(/x). The equations for x|j. on the sub-intervals are coupled due to the terms t — t or 
l-f- 1 - T in (4.2). 



because it will t5Apically have discontinuities at the restarting times t; for i = 1 . . . fc — 1 
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as illustrated in Fig. 4.1 . Thus, we should define the space in which we look for 
solutions as the space oipiecewise continuous functions: 

Q ={x : [—1,0] -^ C" : X continuous on each (half-open) subintervalJ; (4.7) 
and limt/.f x(t) exists for alH = 1 . . . k.} 
Qo={^^Q :x(t,) = Oforalli = 0...fc-l}. (4.8) 

We call X e Q a solution of (4.4)-(4.6) if it satisfies (4.4) for every t e [-1,0]. Both 
spaces, Cj. and Q,o ^r^ equipped with the usual maximum norm. For example, the 
piecewise constant function v^ as defined in (4.5) is an element of Q- Several operations 
are useful when dealing with functions in Q to keep the notation compact. Any element 
X e Q has well-defined one-sided limits at the boundaries t;, which we denote by the 
subscripts -I- and — : 

x(t;)_ = limx(t) fori = l...fc, 

x(t,)+ = lim x(t) = x(t;) for i = . . . fc - 1. 

t\ti 

We also define the four maps 

S[vo...Vfc_i]^(t) = v, if te [t;,t,+i)fori = 0...fc-l, 

r+[x(0] = [x(-l)+,x(tl)+,...,x(t^_l)+]^ 

r_(M)[^(-)] = [MO)_,x(tl)_,...,x(tfc_l)-]^ 
where 

^(s)x(s)+B(s)^^;-:) ^ ';^\'-''\^ 

[nx[l+s -t) ifse[-l,T-l) 

(4.9) 
(the piecewise constant function a^^ was defined in (4.6) as a,t(0 = f; if t ^ [f;, f,+i))- 
The map S takes a tuple of vectors Vo,...V;t-i ^nd maps it to a piecewise constant 
function, assigning x(t,) = V; and then extending with a constant to the subinterval 
J;. For example, the function v^ in (4.4)-(4.5) is equal to Sv. The map F^ takes the 
right-side limits of a piecewise continuous function (thus, r+S is the identity in C"''). 
The map r_(/x) takes the left-side limits at the interior boundaries t; (i > 1), and in its 
first component it takes the left-side limit at t^^ = (the end of the interval) multiplied 
by /I. The map Mj.(|U) is the generalization of Mj defined in (2.9): in the definition of 
Mf. the lower boundary in the integral is not —1 but aj^[t) such that we can estimate its 
norm: 





S 


C" 


^-Q 




r+ 


Q 


^c"*" 


r_ 


-(m) 


Q 


^c"*" 


M 


fc(M) 


Q 


^Ck.o 




[ 


Mkitx)x] (t) 



k 



max IKOII^^-MmI max ||B(t) 

te[-l,0] t€[0,l] 



-!^^pi. (4.10) 



According to the Theorem of Arzela-Ascoli (see [4], page 266) the operator Mj.(/i) is also 
compact because it maps any bounded set into a set with uniformly bounded Lipschitz 
constants. Using the above notation we can write the differential equation (4.2)-(4.3) 
(or, rather, the corresponding integral equation (4.4)) as a fixed point problem in Q. 
For a given v = [vq, . . . , v^^.j] e C"*^ and /i e C we find a function x e Q such that 

x = Sv-l-M,,(|u)x. (4.11) 

Due to the norm estimate (4.10) on M^C/i) we can state the following fact about the 
existence of a fixed point x of (4.11) for a given v and /i. 
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Lemma 4.1 

Let R > 1 be given. If we choose fc > Q(R) then the system (4.2)-(4.3) has a 
unique solution x e Q for all initial values (vg . . . vj._i)^ and all /i satisfying 
M<R. 

Proof:. The differential equation (4.2)-(4.3) is equivalent to fixed point problem 
(4.11). The norm estimate (4.10) implies that / — Mj.(/i) is invertible such that the 
solution X of (4.11) is given by [/ - M^(/i)]"^Sv. D 

Lemma 4.1 implies that for /i satisfying Q(||U|) < fc the operator / — Mj.(|U) is an 
isomorphism on Q. 

4.1. Extended characteristic matrix. We can use Lemma 4.1 to define the ex- 
tended characteristic matrix A^^Cm) for jjl satisfying C^(|/i|) < fc. This is a matrix in 
^[nk)x[nk) ^^^ jg defined as 



Afc(/x)v = A(/i) 



Vjc-l 



Vq - /ix(0) 

_1-X(tfc_l)_ 



(4.12) 



where x = [/ — M^.[fjL)~\ ~ S v e Q is the unique solution of the coupled system of initial- 
value problems, (4.2)-(4.3), with initial condition v = (vg, . . . , v^_i)^ e C"*^. The first 
row in the right side of definition (4.12) corresponds to the boundary condition (2.4). 
The other fc — 1 rows, V; — x(t;)_, guarantee for v e ker A(|a) that the right and the left 
limit of X at the restarting times t, agree, making the gaps shown in Fig. 4.1 zero. Thus, 
X = [/ -M;t(M)]Sv is not only in C^. but in C([-1,0];C") if A^(/i)v = . What is more, 
if A;t(M)v = then the right-hand side of (4.2) for the solution x = [/ - Mj^(^)]~'^Sv 
is continuous such that x is continuously differentiable and satisfies for every t the 
differential equation (2.3)-(2.4), which defines the eigenvalue problem for T. 

Let us list several useful facts about Q, Q q, A^.{fjL), T and the quantities defined in 
(4.9). 

1. We can express A^(/x) using the maps defined in (4.9): 



A;,(M) = /-r_(/i)[/-M,(/x)] 's. 
2. Q is isomorphic to C"*^ x Q o- 



(v, if) e C'"^ X Q 
3. The differential equation 



■ Sv + 1/3 e Q 



xC, 



(4.13) 



(4.14) 



k,0- 



xlit) = Alit)x + Bit) 



x(t-T) ifte[T-l,0], 
^0 if te [-1,T-1), 

x(-l)+ = 
x(t,)+ = x(t,)_ for i = l,...,fc— 1 



(4.15) 



has only the trivial solution (it is identical to (2.5)-(2.6) for fx = 0). Its integral 
formulation is x = Sr_(0)x + M^(0)x. Thus, the operator / - Sr_(0) - M^(0) : Q ^ 
C([-1,0];C") (which is of Fredholm index 0) is invertible. 

4. The time-1 map T of the linear DDE x(t) =A(t)x(t) + B(t)x(t - t), defined 
in (2.2) for initial history segments in x e C([ — 1,0];C"), can be easily extended to 
initial history segments in Q. In fact, the extension of T to Q can be defined by 



T=[i- sr_(o) - Mfc(o)] ' [s(r_(i) - r_(o)) + m^ci) - m,(o)] , 



characteristic matrices for periodic DDEs 9 

which maps Q into C([-1,0];C"), such that 

/ - MT = [/ - Sr_(0) - Mfc(O)] "' [/ - Sr_(M) - M^Cm)] , (4.16) 

which maps Q into Q (see Appendix A for a detailed decomposition of expression 
(4.16)). Since the extension of T maps into C([— 1,0];C") (the original domain of the 
non-extended T) its spectrum is identical to the spectrum of the original T. 

Lemma 4.1 and the above facts permit us to apply the general theory developed in 
[8], relating the extended characteristic matrix A^^Cm) ^ c'^"''-*^'-"'^'' to spectral properties 
of the time-1 map T: eigenvalues, eigenvectors, and the length of their Jordan chains. 

Lemma 4.2 (Jordan chains) 

Assume that < |/ij < R and k > Q(R) where Q(R) is defined by (4.10). If 
A^tC/^t*) is regular then l/fJL.^ is in the resolvent set of the monodromy operator 
T. If A;t(M*) is not regular then A* = l/fJ,^ is an eigenvalue of T. The Jordan 
chain structure for A^ as an eigenvalue of T is also determined by Aj. : 

II. the dimension Zg of kerA^^C/^*) is the geometric multiplicity of A^, 
2. the order p of the pole /x^ of Aj^Cm)"^ is the length of the longest 
Jordan chain associated to A^, 

3. the order of /i^ as a root of det A^(/>() is the algebraic multiplicity of 

^*, 

4. the Jordan chains of length less than or equal to p can be found as 
solutions (jo, ■ ■ . jjp-i) (where yp_i is non-zero) of the linear system 

= Afyo 

= A«yo + Afyi 

0=^yo + A™yi+Afy2 (4-^73 



^(P-i) 
° = (F^^° ^ • ■ ■ ^ ^'^^'-' ^ ^'^^'-' 
where A^^ is the jth derivative of A(.(/i) in /x^ and A^ = A^^C/i^). 
System (4.17) is equivalent to the requirement 

Afc(M) [yo + (/^ - fJ-Jji + . . . + C/i - /^«F" Vp-i] = O ((/i - nJP) 

for all iJLK pL-t (this is, in fact, the general definition of a Jordan chain for a nonlinear 
eigenvalue problem). The vectors jj in the Jordan chains found as solutions of (4.17) 
are still vectors in C"*^. We give the precise relation between the vectors y^ and the 
generalized eigenvectors of / — /iT in Equation (4.21) in the proof below. 

Proof:. The proof extends and modifies the construction used by [13], and makes 
the general theory of [8] applicable to time-periodic delay equations. 

First, we give a brief summary of the relevant statement of [8]. Let G : n -^ 
^{Xii Yi) and H : n -^ -S?(^2; ^2) be two functions which are holomorphic on the open 
subset n of the complex plane C, and map into the space of bounded linear operators 
from one Banach space X^ 2 into another Banach space Yj 2- The functions G and H are 
called equivalent if there exist two functions E -.Q-^ ^{X2;X{) and F : n -^ ^{Ji, ^2) 
(also holomorphic on fi) whose values are isomorphisms such that 

H(/i) = F(/x)G(//)£(/i) for all /x e fi. 
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If two operators G and H are equivalent then every Jordan chain (uo,...,Up_i) of 
the eigenvalue problem H{pL)u = for ju^ corresponds to exactly one Jordan chain 
(vq, . . . , Vp_i) of the eigenvalue problem G[fj,)v = at fx^ via the isomorphism 

Vq = EqUq 

vi = £o"i + EiUq 

(4.18) 

Vp_i = £o"p-i + • • ■ + £p-i"o> 
where the operators Ej are the expansion coefficients of the isomorphism £ in /x = /i*: 

£(/x) = £o + (i" - />f*)£i + . . . + (/x - /xJP-i£p_i + 0((/i - iJL.y). (4.19) 

We construct the equivalence initially only for the case of a single delay t < 1. The 
necessary modifications for arbitrary delays can be found in Appendix B (the underlying 
idea is identical but there is more notational overhead). Also, we restrict ourselves here 
to merely stating what the operators H, F, G and E are, relegating the detailed checks 
and calculations to Appendix A. 

The subset of permissible fx, Q. <z C, is the ball of radius R around 0, where the 
initial-value problem (4.2)-(4.3) has a unique solution. In our equivalence the operator 
functions G and H are 

GiiJ.) = I - iiT : ^i = Q ^ ^1=^1, 



Afc(M) 



H(m) = ' j\: X2 = C"'x Q,o - Y2=X 



2- 



FCm)!/- : 



Note that the spaces Xi and X2 are isomorphic via relation (4.14). The eigenvalue 
problem G(/i)v = is the eigenvalue problem for T (reformulated for /i = A~^). The 
eigenvalue problem H(/x)[v, (/;] = is equivalent to Ai^{^)v = 0, (/? = 0. Thus, establish- 
ing equivalence between G and H proves Lemma 4.2. We define the isomorphisms E 
and F as 

£(m)[v, f]=[l- M,(/i)]-' [Sv + ^] 

>i(m)^" 

^ [r+ - r_(o)] xp + r_(M) [/ - m^^Cm)] "' [i - sr+ - m^w] ip' 
[/-sr+-M^(o)]i/i 

(4.20) 

See (4.13) and (4.16) for expressions giving A^^Cm) and / — fxT, defining G and H in 
terms of S, r± and M^. The inverse of / — M^^C/^) is guaranteed to exist by Lemma 4.1. 
The only choice we make is the definition of E, which is clearly an isomorphism from X2 
to Xj (because (v, ^p) -^ Sv + ^p is an isomorphism from X2 to Xj). The definition of F is 
then uniquely determined if we want to achieve the equivalence F(/i)G(/i)£(|U) = H(/x), 
and follows from a straightforward but technical calculation, given in Appendix A. The 
concrete form of E and F allows us to specify the precise relation between Jordan chains 
of Ai.(n)y = and (/ - fxT)x = 0: (xq, . . .,Xp_i) e C([-1,0];C")p is a Jordan chain 
of (/ — /ir)x = in />( = /i^ if and only if 

^o=-Eo[yo,0], yo=r+Xo 

^1 =-Eo [yi, 0] + £1 [jo, 0] , Ji =r+xi - r+£i [yo, 0] 



(4.21) 



fp-i =£0 [yp-i, 0] -I- ... -I- £p_i [yo, 0], yp_i =T+Xp_i - T+Ej [y^.j, 0] 

-r+£p_i[yo,0]. 
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where (jq, • • • , Jp-i) is a Jordan chain of Aj.(|U) in /x = /^^ and Ej are the expansion 
coefficients of £(/i) in /x^ (see (4.19)). In relation (4.21) we used that the finite- 
dimensional component of £(/i)~^x is F+x (see (A.2) in the Appendix). D 

5. Conclusions. We have generalized the construction proposed by Szalai et al. 
[13] such that it can be used to find a low-dimensional characteristic matrix for any 
Unear DDEs with time-periodic coefficients. 

A shortcoming of our construction is that the function det A^^Cm) does not depend 
smoothly on the delay t. This problem is purely due to the particulars of our construction 
because one expects that the roots of detAj.(-) depend smoothly on t. One can 
remedy this problem by constructing a characteristic matrix A^^CO that depends on the 
Floquet exponents instead of the Floquet multipliers. Then one can operate in the space 
Cpgj.([— 1,0];C") of fc times continuously differentiable functions on [—1,0] that are 
periodic in all derivatives up to order fc. One also has to choose a map S : C"'^ that maps 
into Cpgj.([— 1,0];C"), which requires a modification of the construction of the fixed 
point problem (4.4)-(4.6). 

An open question is how our extended matrix is related to the characteristic matrix 
introduced in the textbook of Hale and Verduyn-Lunel [6] for periodic DDEs with delay 
identical to the period. In particular, what happens to the poles of A(/x) as the delay 
approaches 1 (the period)? 

Appendix A. Details of the equivalence in the proof of Lemma 4.2. 

Inverse of E^fx). Let [v, i/?] = [(vq, .. ., vj._i)-'^, i^] e C"*^ x Q_q be an element of the 
space X2. Then the definition of x = £(/i)[v, (f] in (4.20) means that x is the solution 
of the fixed point problem 

X = Sv + (f + Mf.[^.)x . (A.l) 

Consequently, we can recover v by applying r+ to (A.l) (F+i/j = since ^p e Qg ^^d 
r+Mjf(/x) = since M^. maps into Q,o)- Then i^ can be recovered as y) = x — SF+x — 
M^(/i)x such that the inverse of £ is 



£"^(m)x : 



F+x 
[/-SF+-M,(/x)]x 



(A.2) 



Expression (4.16) /or / — /xT. Before finding an expression for F(/x) we check that 
G(|U) = / — pLT is indeed given by the expression (4.16). If we denote the image of x 
under G(/x) by y then y has the form x — y where y = piTx. By definition (2.2) of 
the monodromy operator T, if x e C([— 1, 0]; C") then y is the unique solution of the 
inhomogeneous initial- value problem on [—1,0] 

Ko=A(,»(.)+B(op''-:' , ';;!J7''°' (a.3) 

[[ix{l + t — T) ifte[— 1,T — Ij 
y(-l) = M^(0) 
Kt,)+ = Kf,)-fori = l...fc-l. 

This solution y is also in C([— 1,0];C"). We note that the initial conditions in (A. 4) 
mean that F+y = F_(0)y + [F_(/i) — F_(0)]x. Thus, the fixed point formulation 
corresponding to (A.3)-(A.4) is 

y = S [F_(0)y + [F_(m) - F_(0)]x] + M,(0)y + [M,(m) - M;,(0)]x, (A.5) 

which can be rearranged for y because the factor / — SF_(0) — Mj^^O) (which is in front 
of y) is invertible as explained below (4.15) in the hst of useful facts in §4: 

y=[l- SF_(0) - Mfc(O)] -' [SF_(m) - SF_(0) + M,(/i) - M,(0)] x. (A.6) 
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We notice that the operator in front of x is also defined for all x e Q, which allows us 
to extend the definition of x -^ y = fiTx to x e Q. Moreover, when applying r+ to 
(A.5) we observe that y is continuous in [—1,0], that is, the map x -^ fiTx maps Q 
into C([— 1,0];C"). Inserting y = x — y into (A. 6) gives expression (4.16), which is 
equivalent to the fixed point problem 

y = Sr_(0)y + M^Wy + x- Sr_[^)x - M,(/i)x. (A. 7) 

Determination ofF^fx). From the fixed point equation (A. 7) it becomes clear how 
to choose the other isomorphism F(/x) such that F(/i)G(/x)£(/x)[v, c^] = H(/i)[v, if]. We 
observe that (A. 7) implies (by applying r+ to both sides) 

r+y = r_(o)y + r+x-r_(M)x. (a.s) 

We apply S to this identity and subtract it from the fixed point equation (A. 7) defining 

y- 

y - Sr+y = M^(_0)y +x- Sr+x - M^(m)x. 

If X = £(/x)[v, (/)] = [/ -M;,(/i)]"^[Sv + t^] then (A.2) implies T+x = v, such that we 
obtain for y = G(/>i)x 

y-Sr+y-Mfc(0)y = (^. (A.9) 

Thus, [/ -Sr+ - M^(0)]G(/x)£(/i)[v,(/7] = if for all (v, (/?) e X^, which justifies our 
choice of the second component of F(/x) in (4.20). Inserting the expression £(/x)[v, if] = 
[I - Mi^(ifi)]~^ [Sv + if] for X and r+x = v into (A.S) we obtain that 

r+y - r_(0)y = V - r_(/i)[/ - M;,(m)]-' [Sv + if] (A.10) 

= [/ - r_(M)[/ - m,(m)]-is]v - r_(M)[/ - m,(m)]-V 

= A^(/i)v - r_(/i)[/ - Mfc(M)]-V, (A.11) 

where we used the definition (4.13) of the characteristic matrix A^^Cm)- Since we know 
that If can be recovered from y via (A.9) we obtain from (A. 11) the first component of 
the definition of F(|U) in (4.20): 

r+y - r_(0)y + r_(M)[/ - M.C/i)]-^ [/ - Sr+ - M,(0)] y = A^[fi)v. 

Inverse of F(/i). It remains to be checked that F is invertible if |/x| < R. Given 
V e C"*" and if e Q q such that 

V = r+y - r_(o)y + r_(M)[/ - m,(/x)]- V, (A.12) 

i^ = y-Sr+y-Mfc(0)y, (A.13) 

how do we recover y? Applying S to (A.12) and adding the result to (A.13) gives 

Sv + if= Sr_(M)[/ - M;t(M)]" V +[I- Sr_(0) - M^CO)] y, (A.14) 

which is equivalent to (A.12)-(A.13) because the map [v, f] e C"'' x Q,o -^ Sv + (/? e Q 
is an isomorphism. Equation (A.14) can be rearranged for y because / — Sr_(0) — Mj.(0) 
is invertible. Hence, [v,if] e C"*^ x Qq and y e Q satisfy [v, (/;] = F(/i)y (that is, 
system (A.12)-(A.13)) if and only if 

y = [/ - Sr_(0) - M,(0)] -' [Sv + f- Sr_[fi)[I - Mfc(/x)]- V] ■ 
Consequently, 

f (m)"' [v, (^] = [/ - sr_(o) - M,(o)] -' [sv + v^ - sr_(/i)[/ - M,(/i)]- V] ■ 

(A.15) 
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Appendix B. Characteristic matrices for general periodic delay equations. Con- 
sider a linear DDE with coefficients of time period 1 and with maximal delay less than 
or equal to an integer m > 1: 



x(t)= d0r?(t,0)x(t-0) (B.l) 

10 

where r]{t,6) is bounded, measurable and periodic in its first argument t (with period 1), 
and of bounded variation in its second argument 6. We choose rjCt, •) e JVBV(K; C"^"), 
that is, r]{t, •) = on (— oo, 0], r]{t, •) is continuous from the right on (0, m), and 'q^t, ■) 
is constant on [m, oo). In addition we assume that the total variations of all Tj(t, •) have 
a common upper bound: 

V^"rj(t,0<V (B.2) 

where V^f is the total variation of/ [2]. In formulation (B.l) the dependent variable 
is the function x on the history interval [— m, 0]. The time-1 map Tx of (B.l) for 
X e C([-m, 0];C") is defined as the solution y e C([-m, 0];C") of the equation 



yit-e) ift-0>-l 

t<-l: y(t) = x(l+t). (B.4) 

Only the part of y on the interval (—1,0] is the result of a differential equation (namely 
(B.3)). The remainder of y, y|[_„_i], is merely a shift of x, given by (B.4). Note 
that we have included the initial condition for (B.3) into (B.4), namely y(— 1) = x(0). 
Correspondingly, the eigenvalue problem for T reads 



. ^ r . ^^ (^(t-0) ift-0>-l 

'>-'■■ ^^^) = J/^^^^'^nMx(l + t-.) ifr-.;-l '^-'^ 

t<-l: x(t) = /ix(l+ 1). (B.6) 

The spaces Q and Q Q' ^^d the maps S, Fj- and M^.[^) defined in (4.7), (4.8) and (4.9) 
for a single delay can be extended in a straightforward way. We partition the interval 
[-m,0] into sub-intervals of length 1/fc, J, = [t,,t;+i) = [-1 + i/fc,-l + (i + l)/fc) 
(i = fc — mk ...k—1; note that i can become negative), and define the spaces 

Q = {x : [— m, 0] -^ C" : X continuous on all (half-open) subintervals J,, 

and limf^t x(t) exists for all i = fc — mk + 1 . . . fc} (B.7) 
Q ={x e Q : x(t,) = for all i = . . . fc - 1}. 

Note that non-negative indices i of t; correspond to times t; in the interval [—1,0] 
whereas negative indices i correspond to times t; in the interval [— m, — 1). These points 
of discontinuity are treated differently in the definition of Qq- to be an element of 
Q,o the function x is only required to be zero at t; with non-negative i. The space Q 
consists of piecewise continuous functions where the points of discontinuity are the 
times t, e [— m, 0]. For functions in Q we define the maps (again analogous to the 
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definitions in (4.9) for the single delay) 



^nk 



S[Vo...Vt_i]^(t): 



r+ 


Q -^ C"'' 


r_(/x) 


Q -^ C"'' 


M;t(M) 


Q '^ Q,o 



v; if te [t;,t,+i)fori: 
ifte[-m,-l), 

r+x = [x(-l)+,x(ti)+, . . .,x(t^_i)+] ^ , 
r_(/i)x = [Mx(0)_,x(ti)_, . . .,x(tfc_i)_] ^ , 
where 
( 



0...fc-l, 



[M;t(M)x] (0 ■■ 



\ 



J JdeT,(s,0) 



ak(t) 



x(s - 0) 



/ix(l +s — 0) 



if s — i 
if s — i 



Mx(l + t) 



>-l 
<-l 

for t > -1, 
for t < -1, 

(B.8) 
and as.(t) = t; if t e J; and i > 0. The map S extends a tuple of vectors into a piecewise 
constant function using the elements of the tuple as the values on the subintervals 
J; c [—1,0], and setting the function to on [— m, — 1). T^ and r_ are defined in 
exactly the same way as in (4.9) (again, the notation x(t;)± refers to left-sided and 
right-sided limits of x at t;). The definition of M^^Cm) is identical to the single delay case 
for t > —1. For t < — 1 we define M^(/x) as a shift multiplied by /i. In the same manner 
as for the single delay we consider the extended initial-value problem 



t >-l : 



t<-l : 



X(t): 



deVit,e) 



x(t-0) 



if t-0>-l, 



(B.9) 



[nx[l + t-9) ift-0<-l, 
x(t;)+ = V; fori = 0...fc-l, (B.IO) 

x(t) = Mx(l-|-t). (B.ll) 



Note that the initial values v, are given only at the discontinuity points t, e [—1,0] 
(i > 0). The solution x e [— m, — 1) is defined by a backward shift using (B.ll). The 
precise meaning of the initial-value problem is given by its corresponding integral 
equation, which is a fixed point problem in C^^, and can be expressed as 



X =Sv -|-Mj.(/i,)x 



(B.12) 



using the maps S and M^^C/x). The following lemma gives a sufficient condition for the 
existence of a unique solution x e Q of (B.12) and the invertibility of / — Mj.(/x). It 
is a straightforward generalization of Lemma 4.1. The only difference in the proof is 
that we can achieve ||Ms.(/i)|| < 1 only in a weighted maximum norm || • ||^ (which is 
nevertheless equivalent to || • H^o). 

Lemma B.1 (Unique solution for initial-value problem) 

LetR > 1 be fixed. Then the fixed point problem (B.12) has a unique solution x 
for all V e C"*^ and for all complex /i satisfying \pL\ < R if the number k of the 
sub-intervals is chosen such that 



pm+1 '-' 



(B.13) 



Note that V is an upper bound on the total variation of 17. The left-hand side of (B.13) 
approaches 1 for fc ^ 00 whereas the right-hand side is larger than 1 if R > 1 and 
R> |/i|. 

Proof:. We define the norm 

||x||jj= max |R^x(t)| 

te[-m,0] 
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on the space L°°([— m,0]; C") of bounded functions and the spaces Q and Q q (which 
are subspaces of L°°). This norm is equivalent to the standard maximum norm. Thus, 
(B.12) has a unique solution in Q whenever ||M;t(/^)llR < 1- Two operators that appear 
in Mj^^fi) are shift and integration. Their respective domains and norms are: 



S(0):L°°([-m,O];C")-^L°°([-m,O];C") 

x(t-0) if t-0e[-m,O], 
otherwise, 

J:L~([-m,0];C")-Qn 



[S(0)x](t) 



[Jx](t)= x(s)ds. 



defined by 
||S(0)x||«<R'^||x||„ 
defined by 

||Jx||„< ||x|| 

" "^- logR " " 



In order to find the || -Hjj-norm of the integrand (which is in L°°([— m, 0] ; C") if we extend 
it with zero for t e [— m, —1)) in the definition of M^.{pL), (B.8), we have to estimate the 
II ■ ||j;-norm for a Stieltjes sum over an arbitrary partition < 6q < . . . < 6;^ < m (with 
arbitrary intermediate values Oj e [9j, 0j+i]): 



N-l 



sr^r . ^ . . ^ol^(f-0,) ift-0,>-l 



N-l 

< sup R'yUt,ej+,)-r]it,ej)\ 

t6[-l,0] p^ 

N-l 



[S[9j)x] (t) 
[S(0, - l)x] (t) 



if t-0j >-l 
if t-0, <-l 



t€[-mj^' ' I (R||S(0j-1)x||r ift-0j<-l 



l|S(0j)x||jj 



if t-0j >-l 



JV-l 



< sup V|T}(t,0^+i)-r,(t,0j)|R^^||x||« 
te[-l,0] j=o 

<R"'\\xU sup V-r,(t,0 = il'"||x||«V. 

t€[-l,0] 

Combined with the norm estimate for integration from a^.[t) to t we obtain that 



[M;,(m)x] 



[-1,0] IIr 



< 



lUR-VllxlUSR-V- 



.Ri/''- 



logR 



The part of M^(/i)x on the interval [— m, — 1) has the || • H^-norm: 

II [M,(/i)x] |[_,„._i)||fi < |m|||S(-1)x||« < \ii\R-'\\xU, 



(B.14) 



(B.15) 



where we note that S(— 1) maps Q back into Q- Adding up the inequalities (B.14) and 
(B.15) we obtain an upper bound for ||M;^(|U)x||s: 



l|Mfc(M)x|U < 



R'^V + — 

logR R 



(B.16) 



The factor in front of ||x||s on the right-hand side is less than unity if k satisfies the 
inequality (B.13) required in the statement of the lemma. D 

The invertibility of / — M^(/x) permits us to choose exactly the same constructions for 
the characteristic matrix A^^Cm) ^nd for the isomorphisms E and F proving Lemma 4.2 
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in the same way as for the single delay case (see (4.13), (4.16) and (4.20)): 

G(m) = I-fjiT=[l- Sr_(0) - Mfc(O)] "' [/ - Sr_(M)Mfc(/i)] , 

A;,(/i) = / - r_(|a) [/ - M^(/i)] "^ S e C"'' x C''^ 
£(/x)=[/-M,(/i)]-'[Sv + (^], 

"r+ - r_(o) + r_(/i) [/ - m^(m)] "' [/ - sr+ - m^co)] ,' 



f(m): 



/-sr+-M,(o). 



With these constructions the maps £(/x) : X2 -^ X^ and F(/>() : Xj -^ Xj are isomorphisms 
for |/i| < R (see (A.2) and (A.15) for the inverses). The relation H^fi) = F(/i)G(/i)£(/i) 
makes the infinite-dimensional eigenvalue problem of the time-1 map, [/ — fj,T]x = 0, 
equivalent to the finite-dimensional eigenvalue problem Aj.(/^)v = 0. 

We remark that the inverse of / — Sr_(0) — M^(0), which is present in the definition 
of G(/i), exists because the initial-value problem (B.3)-(B.4) has only the trivial solution 
for X = (the operator is a compact perturbation of the identity). In other words, the 
existence of this inverse follows from the fact that the time-1 map Tx is well defined 
(and equal to zero) in x = 0. Moreover, the extension of T to Q, defined by 

T=[i- sr_(o) - Mfc(o)] "' [s(r_(i) - r_(o)) + m^(i) - m^(o)] , 

which maps Q back into Q, has the property that T^x is continuous (that is, r'"x e 
Ci[—m, 0]; C")) for all x e Q- Hence, the spectrum of the extension T : Q '^ Q is the 
same as the spectrum of the original T : C([-m, 0];C") -^ C([-m, 0];C"). 
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